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ABSTRACT 


A  survey  is  made  of  the  mathematical  properties  of,  and  the 
arithmetic  relationships  between  various  distributions  on  the 
circle  and  the  sphere.  Some  new  results  are  given.  The  Brownian 
motion  and  Angular  Gaussian  distributions  are  shown  In  computer 
drawn  graphs  to  bracket  the  von  Mises-Fisher  distributions. 
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1.  Introduction 


Because  I  began  my  work  on  directional  data  while  in 
Professor  Moran's  Department,  It  seems  appropriate  to  take 
up  the  topic  again  for  this  volume  in  his  honour.  My  motivation 
then  was  palaeomagnetism  which  has  since  greatly  grown  in 
scientific  importance.  The  original  direction  of  magnetization 
of  rocks  proved  to  be  an  Invaluable  clue  to  the  discovery 
of  plate  tectonics.  Later,  discussions  with  ornithologists 
about  the  navigational  and  homing  ability  of  birds  led  me 
to  an  Interest  in  directions  in  two  dimensions.  In  the  inter¬ 
vening  25  years,  many  statisticians  have  created  a  sizeable 
literature,  and  a  book  (Mardta,  1972),  on  the  statistics 
of  directions.  Though  the  theory  often  applies  to  an  arbitrary 
number  of  dimensions,  we  will  consider  only  2  and  3  dimensions. 

Almost  all  the  parametric  work  has  been  based  on  one 
distributional  type,  introduced  in  another  context  by  von  Mises 
(1918)  for  the  circle  and  developed  for  palaeomagnetism  by 


Fisher  (1953). 


The  nomenclature  is  unfortunate  because  Arnold  (1941)  in 
an  unpublished  M.I.T.  Ph.D.  Thesis  derived  the  distribution 
named  after  Fisher  and  gave  the  maximum  likelihood  estimation 
procedures  for  both  the  von  Mises  and  Fisher  distributions. 

He  also  derived  the  Brownian  motion  densities  (6)  and  (7). 

The  von  Mises  density  is 


,m(e-K>’ exp (iccose) 

tfiere  tc  a  o  and  IQ(<)  Is  a  Bessel  function, 
direction  r  makes  an  angle  0  with  the  "mean" 


(1) 

The  observed 
direction 
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The  Fisher  density  is,  with  k>0, 

*  4 tt s i n H~ exp  (<cos0)  (2) 

using  spherical  polar  coordinates  (6,<|>)  so  that  r*y=cos0.  In  both 
cases,  k  is  a  "concentration"  parameter.  When  tc  *  o,  the  density 
is  constant.  When  <  is  large,  the  probability  is  concentrated 
about  0*0.  The  exponential  form  makes  inference  procedures 

easy  to  develop.  It  is  symmetric  about  the  mean,  or  polar, 
direction.  It  shares  some  statistical  characterizations  of  the 
Gaussian  distribution,  which  it  is,  approximately,  for  large  k. 

See  e.g.  Mardia  (1975).  More  important, the  exponential  form  of 
(1)  and  (2)  makes  statistical  inference  procedures  fairly  easy  to 
develop. 

Multipolar  forms  may  be  created  by  replacing  0  and  by  p0, 
p  an  integer.  When  p*2,  the  density  is  proportional  to 

p 

exp  (k  cos  0)  (3) 

and  gives  bipolar  (<  >  o)  and  girdle  distributions  (k  <  o)  on 
the  sphere.  Its  generalization  has  density  proportional  to 

exp  (kj  Xj2  +  k 2  x22  +  k3  x32)  (4) 

where  x  *  (xj,  x2,  x3)  Is  the  observed  unit  vector  referred  to 

three  orthogonal  axes.  If  the  data  Is  x^,...x^,  then  the 

statistics  of  (3)  and  (4)  depend  upon  the  eigen  values  and  vectors 

of  Z  ,(Ox(ir.  whereas  the  statistics  of  (1)  and  (2)  are  based  on 
1 

n  Ml 

I  xv  Here  x  Is  a  column  vector,  %'  Its  transpose. 

1 

Fisher  pointed  out  that  the  distributions  (1)  and  (2)  can  be 


3 


obtained  from  a  multivariate  Gaussian  by  conditioning  as  will 

be  explained  In  the  beginning  of  52.  The  distributions  (3)  and 
(4)  also  come  from  a  Fisherian  construction.  Let  XltX2>X3  have 
the  joint  density 

(2ir) ’3/2(a1  a2  ajJ’^exp-l^x^/oJ+Xg/o^Xj/o*) 

If  R2  ■  X2  +  X^  +  X2  then  the  distribution  of  the  direction  of 

X  when  R  Is  conditioned  to  be  unity  Is  the  Bingham  distribu¬ 
tion.  If  Oj  *  a 2 ,  this  reduces  to  (3).  If  °i  “  °2  >  >  °3» 
clearly  makes  (3)  a  girdle  distribution  with  the  probability  con¬ 
centrated  around  a  great  circle  In  the  x^-x^  plane. 

There  has  always  been  the  possibility  of  constructing  den¬ 
sities  which  are  the  exponentials  of  higher  order  polynomials  In 
the  components  of  x.  It  Is  natural  to  rearrange  these  to  obtain 

exp  ZBpSp(0 ,♦) 

where  $m(e,$)  are  surface  harmonics  and  8  &  ♦  are  spherical  polar 
angles.  Recently  Beran  (1979)  has  suggested  how,  despite  the 
complicated  normlng  factor,  to  find  maximum  likelihood  estimates 
of  the  ®p's.  Our  efforts  to  apply  his  method  will  be  reported  in 
Watson  (1981).  For  p>2,  these  distributions  cannot  have  a  Gaussian 
origin.  It  Is  not  known  whether  they  could  arise  from  diffusions 
of  the  types  mentioned  below. 

On  the  unit  circle,  one  may  get  distributions  by  "rolling 
up"  distributions  on  the  line.  This  was  a  topic  of  Interest  to 
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mathematicians  In  the  30's  and  40's  -  see  e.g.  Levy  (1939), 

Havlland  (1941),  Wlntner  (1947).  Arnold  (1941)  In  his  thesis 

considered  this  process  and  diffusion  and  suggested  (1),  (2),  (3) 

and  (6).  Thus  a  density  f(x)  on  the  line  becomes,  on  the  unit 
circle. 


h(x)  f  (x  +  2-n p ) . 

If  the  Gaussian  Is  used,  this  procedure  leads  to  the  Brownian 
density  (6)  mentioned  in  the  next  paragraph.  The  procedure  has 
no  analogue  on  the  sphere  since  It  Is  not  a  developable  surface.  It  Is 
an  Interesting  conjecture  that  distributions  on  ruled  surfaces 
will,  one  day,  be  generated  from  linear  densities. 

Another  simple  process  for  generating  circular  from  linear 
or  planar  distributions  Is  to  use  bl-llnear  (or  linear  fractional) 
transformations.  This  suggests  e.g.  putting  x  »  tan  e/2  which 
sends  the  standard  cauchy  Into  the  uniform  distribution.  There 
are  also  some  connections  with  number  theory  -  see  e.g.  Stapleton 
(1963)  and  Schmidt  (1977).  We  will  discuss  these  on  another 
occasion.  We  turn  now  to  the  major  other  source  of  distributions  - 
diffusion  processes. 

If  a  particle  starts  at  8  ■  0  and  executes  symmetrical 

Brownian  motion  on  the  circumference  of  a  unit  circle  with  a 
2 

variance  o  per  unit  time,  the  density  at  e  at  time  t  Is 

_  1  f  i  a  O _  \2  1 


br2(e,v)  »  ?  —J ==—  exp  j- 

c  p.  .os  o/2irt  2o  t  j 

•  ^1  ♦  2  ?  exp  (-£^v)cos{  p  0  )j 


(6) 
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2 

where  V  ■  at.  On  the  line  the  density  would  be  Gaussian  with 
mean  zero  and  variance  V;  this  density  when  substituted  In  (5) 
gives  an  alternative  form  of  (6).  A  symmetrical  Brownian  motion, 
starting  at  the  6-0  point  of  a  spherical  coordinate  system, 
yields  a  density  on  a  sphere  of  unit  radius  given  by 


br3(0,$,v) 


1 

Ttt 


? 

p-0 


(2p  +  1)  exp( -p(p+l )V/4) 


Pp(cos  e) 


(7) 


where  Pp(cos  0)  Is  the  Pth  Legendre  polynomial  In  cos  0.  These 

distributions  are  the  probabilistic,  rather  than  statistical, 
Gaussian  analogues  of  the  circle  and  sphere.  In  a  later  section 
we  will  see  other  diffusion  processes  which  lead  to  the  dls- 
trlbutlons  mentioned  earlier. 

Some  other  distributions  have  been  suggested. 

Scattering  theory  of  course  provides  them  In  physics.  Kendall 
(1974)  studied  a  model  of  bird  navigation  and  derived  the 
(complicated)  distribution  of  the  angle  where  the  bird  first 
hits  Its  (circular)  home  territory  -  the  Lack  distribution. 

In  passing  he  discusses  briefly  the  off-set  or  displaced  normal 
distribution  which  has  arisen  in  other  contexts.  We  will  call  It 
the  angular  Gaussian  In  the  next  section.  The  angular  Gaussian 
has  been  suggested  In  two  studies  of  the  precision  of  palaeo- 
magnetlc  measurements  --  Harrison  (1980),  Brlden  and  Arthur  (1980), 
the  latter  paper  providing  approximate  statistical,  methods. 


The  mathematical  and,  particularly,  the  arithmetic  re¬ 
lationships  between  the  von  Mlses-Flsher ,  Brownian  motion  and 
angular  Gaussian  distributions  Is  the  subject  of  this  paper. 

2.  The  angular  Gaussian 

Fisher  observed  that  his  distribution  may  be  derived  from 

a  Gaussian  by  conditioning.  If  X  has  a  trivariate  Gaussian 

2 

distribution  with  mean  vector  y  and  covariance  matrix  o  I,  set 
R  *  | X |  ,  L  *  X/R,  A«y/|y| . 

Then  the  joint  density  of  R  and  L  Is 

(2 


exp(-  A,  ,«*♦  2>}.xp 


so  the  density  of  L  on  the  unit  sphere,  conditional  upon  fixed 
R  has  the  form  (2)  with  <  *  R|y|/o  ,  cos  0  »  t '  X.  It  Is  more 
surprising  that  he  did  not  comment  on  the  way  palaeomagnetlc 
directions  are  actually  obtained.  The  three  components  of  the 
magnetic  Intensity  (let  It  be  X)  are  measured  by  rotating  the 
specimen  In  the  same  apparatus.  Hence  the  above  distribution 
for  X  Is  appropriate  but  the  unconditional  distribution  of 
L  Is  what  Is  required.'  Setting  S  «  R/o  and  m  *  |y|/o,  this  Is 


#93  (L,X,m)  ■ 


S  exp  < 


■  7  (s 


2+n2-2S«L 


'l 


dS  (8) 


It  Is  a  function  only  of  l1^  ■  cos  6  and  so  rotatlonally 
symmetric  about  the  mean  direction  m  Is  an  a  concentration  para¬ 
meter.  The  statistical  aspects  of  (8),  and  Its  relation  to 
palaeomagnetlc  work,  are  discussed  In  Watson  (1980).  Of  course 


If  measurement  error  Is  not  the  only  source  of  variability  in  the 
components  of  X  ,  X  may  well  have  a  general  Gaussian  distribu¬ 
tion  and  distributions  more  complicated  than  (8)  and  (9)  arise. 

00 

ag2(L,X,m)  *  ^  S  exp  jj  j  (S2  +  m2  -2SmL.'X)j  ds  (9) 

The  densities  (8)  and  (9)  give  the  angular  Gaussian  distributions 
on  the  sphere  and  circle  respectively. 

The  Fourier  expansion  for  the  von  Mises  density  (1)  is 

vm{e,K)«  y-  J.  T^-k)  exp  ( 1  p 6 )  (10) 

I  (  1C  J 

“  iir (1  +  2 1  i^r<)  cos  (p0> 

since  Ipl«)  *  I_p(<).  Thus  the  Fourier  expansion  of  ag2(L,X,m) 
with  cos  0  ■  L'x  Is  therefore 

00 

ag2(L,A,m)  ■  e  ^  |  exp  j"  J  ^  +  m  )►  SIp(Sm)dS, 

or  r  0  -v  - 

ag2(e,m)  ■  <  ^(m)  +  2  ?  g^  ^(m)  cos(p0)>  ,  (11) 


where 


•  «f,(.)  ■  s  lp(s"0  exp  <  -  J  (S2  +  m2)  1 

i  irrsmy  L  J 


0  *0 


Similarly,  the  spherical  harmonic  expansion  of  ag^U.x.m)  or 
ag3(e.$.m)  follows  from  the  expansion  of  (2).  This  like  (7)  was 
given  by  Roberts  and  Ursell  (1960)  and  is 


8 


■  ?  2p  +  1  I  +  ,  ( k)  P  (cose).  (13) 

0  ^r~ 

Hence  the  expansion  of  (8)  Is 

a*,e, ♦,««)-  ?  2p  +  1  g£3)(m)  P_(cose)  (14) 

4  0  4tt  p  p 


where 


00 

r 


S2exp  -  y(S2  +  m2) 


0 


sinhSm 

5m 


dS 


The  angular  Gaussian  Is  a  special  case  of  distributions 
obtained  by  integrating  the  length  of  a  vector.  Thus  let  f(z) 
be  the  density  function  of  vectors  z  in  Rq  and  set  z  *  y  l 
where  y  Is  the  length  of  z.  Then  the  distribution  of  the  direction 
l  on  «p,  the  surface  of  the  unit  sphere  is  Rq,  is 


9(1) 


.CO 

yq_1  f(yl)  dy 
o 


(15) 


This  family  of  distributions  is  examined  in  Watson  (1981). 
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3.  Some  mathematical  relations  between  distributions 

The  classical  way  of  matching  two  distributions  with  the  same 
mean  on  the  line  is  to  adjust  the  parameters  so  they  also  have  the 
same  variance.  It  is  easy  mathematically  but  need  not  lead  to  a 
good  probabilistic  match.  For  directions  the  analogue  is  to  match 
E  cos*  8  i.e.  make  the  p  ■  1  terms  in  the  Fourier  (or  spherical) 
expansions  equal.  Now  that  computing  is  easy,  this  or  other 
methods  may  be  used  and  compared.  As  we  will  see,  it  is  not 
possible  to  Improve  much  on  E  cos  0  matching  for  several  defini¬ 
tions  of  a  good  match. 

Roberts  &  Ursell  (1960)  began  this  work  by  showing  that 
the  Brownian  distribution  (7)  and  the  Fisher  distribution  (2) 
matched  well  If  one  chose  <  and  V  so  that  E  cos  0  was  equal 
for  both.  Comparing  (7)  and  (13)  this  means  they  are  related  by 


-  *  exp  (-V/2)  (16 

l/2(< 


Since  both  are  rotationally  symmetric,  the  merit  of  the  matching 
could  be  judged  by  s up (  FI(0,tc)-BR3(0,V)  |  where  F I ( 6 , < )  and  BR3(0,V) 

are  cumulative  distributions  of  0  easily  obtained  from  (2)  and 
(7)  respectively.  They  found  that  then  the  supremum  never 
exceeds  (slightly  over)  .02  which  it  takes  when  exp-V/2  «  1/2. 

The  author  suggested  to  Stephens  that  he  do  the  same  in 
his  thesis  for  the  circular  distributions  (1)  and  (6).  The 
matching  relation  is ,  comparing  (6)  and  (10), 


10 


r^j  *  exp  <-v/2>  07) 

Stephens  (1963)  found  that  the  sup  |VM(8,ic)  -  B Rg ( 0 • V ) I  never 
exceeded  .0125  which  was  attained  for  the  matched  pair  <  =  1.4, 

V  *  1.09.  Further,  he  found  that  the  matching  could  hardly  be 
improved  beyond  that  implied  by  (17).  The  relationships  (16) 
and  (17)  may  be  simplified  when  <  o  and  <  00  . 

These  fascinating  results  led  to  a  number  of  papers.  Hart¬ 
man  and  Matson  (1974)  showed  that  there  is,  in  any  number  of 
dimensions,  a  distribution  depending  on  <  and  the  dimensionality, 
for  the  diffusion  time  t  so  that  the  time  mixture  of  the  Brownian 
distribution  is  exactly  equal  to  the  von  Mises-Fisher  distribu¬ 
tion.  The  Brownian  distributions  are  infinitely  divisible.  Kent 
(1977,  1979,  1980)  showed  that  this  mixing  distribution  is  in¬ 
finitely  divisible  so  that  this  is  also  true  of  the  von  Mises- 
Fisher  distribution.  J.  Pitman  (personal  communication)  has  shewn 
that  it  is  not  unique. 

Kent  (1978  )  showed  that  the  von  Mises(Fisher)  distribution 
is  the  equilibrium  distribution  of  the  analogue  of  the  Ornstein- 
Uhlenbeck  process  on  the  circle  (sphere).  The  generator  of  this 
diffusion  on  the  circle,  e.g.,  is 

G  =  \  dW  -  *  sin  9  d/de  .  (18) 

Reuter  (unpublished,  Kent  1978)  showed  that  the  von  Mises(Flsher) 
arises  if  a  diffusion  with  drift  starts  at  the  origin  and  stops 
when  it  hits  the  unit  circle  (sphere). 

Recently,  Yor  (1979)  has  given  the  explicit  formula  for  a 
mixing  distribution  and  a  probabilistic  explanation  using  two 


-11- 


dimensional  Brownian  motion.  The  formula  was  first  obtained  by 
Kendall  and  appears  in  Kent's  Ph.D.  Thesis  (1976)  but  was  not 
published.  From  (6)  and  (10),  the  von  Mises  will  be  a  mixture 
of  Brownians  if  and  only  if  for  all  integers  p 


y«> 

i0(<) 


■ ./ 


-  e-p2V2 


dnK(V) 


(19) 


where  n^V)  is  a  probability  distribution.  Yor  and  Kendall  showed 


that 


dn  (V)  =  3/  Js 

*  it  /  2  (2V)* 


exp(Tr2/2V)  \|>  (V)dV 


(20) 


I0(<) 


where 


♦«(»>  ■ 


fco 

expi 

o 


■]} 


(-x2/2V)-KCOshx| \s1nhx  simrxdx.  dx.  (20') 


Hartman  and  Watson  (ibid)  give  a  formula  for  obtaining  from 
(20)  a  mixing  distribution  in  any  .(umber  of  dimensions. 

The  angular  Gaussian  is,  like  the  von  Mises-Fi sher ,  a  mixture 
of  Brownians.  For  the  circle,  the  density  is 

1  l  CO  ^ 

ag2  (6,m)  *  ZtT  ^1+2  £  cp(m)cosp6j  , 


where 


exp  £j(s2+m2))  Ip(sn)sds 


(21) 


(22) 


Using  the  Kendall-Yor  formula  for  Ip(sm)  and  Inverting  the 
order  of  integration,  we  see  Immediately  that  the  mixing  density 
for  the  angular  Gaussian  Is 
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I 


oo 

0 


exp  j(s2+m2)j 


^SB(ulIo(s")sds 

du 


(23) 


It  i^  a  density  because 

exp  L  |-(s2+m2)J  s  IQ(ms)  (24) 

is  a  density  —  for  the  length  of  a  normal  vector,  in 

fact.  Similar  results  are  easily  obtained  for  ag^. 

Further,  the  angular  Gaussian  distributions  may  be  infinitely 

divisible,  because  the  further  mixing  distributions  (24)  and 

2 

its  higher  dimension  analogues  are  infinitely  divisible.  R  is 
infinitely  divisible  but  we  have  here  mixed  with  R! 

Returning  to  (18),  we  see  that  Kent's  construction  is  of 
a  circular  variant  of  the  Orstei n-Uhl enbeck  process  -  the  restoring 
force  towards  the  origin  is  proportional  to  sine,  and  his 
equilibrium  distribution  is  simply  the  Boltzman  distribution  since  the 
potential  energy  is  proportional  to  cos9.  Formally,  the  equi¬ 
librium  density  then  satisfies 

\  ^7  "  W~  (Xsin0f)  "  0  (25) 

so  that  the  density  f  is  proportional  to 
exp  (2Xcos8  ) 


If,  however,  the  restoring  force  is  proportional  to  sin26  ,  we  get 
attraction  (X  >  o)  to  9  =  it  as  well  as  9  =  0,  and  the  equilibrium 
density  is  proportional  to 
exp(Xcos29)  , 
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which  is  the  circular  variant  of  (3).  A  more  elaborate 
construction  gives  the  Bingham  distribution  (4). 


d 
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4.  Arithmetic  and  Graphical  Comparisons  of  Distributions 

The  definitions  of  the  von  Mises-Fisher,  Brownian  and  angular 
Gaussian  densities  only  suggest  that  they  may  be  numerically  similar 
when  their  concentration  parameters  are  either  very  small  or  very  large. 
If  they  can  always  be  matched  very  closely  by  appropriate  choices 
of  their  concentration  parameters,  then  for  statistical  analysis  it 
should  matter  little  which  is  used  even  if  a  specific  family  is  in¬ 
dicated  by  stochastic  modelling.  We  may  then  use  the  distribution 
which  is  most  convenient.  For  example  it  is  much  easier  to  simulate 
the  Brownian  and  angular  Gaussian  on  the  circle  than  the  von  Mises 
distribution  since  they  require  only  the  generation  of  one  and  two 
Gaussians  respectively  and  simple  algorithms.  On  the  sphere,  the  angular 
Gaussian  and  Fisher  distributions  are  easier  than  the  Brownian  since 
they  require  only  three  Gaussians  and  three  uniforms,  respectively,  an«* 
a  simple  algorithms.  For  statistical  mathematics  other  "rankings" 
obtain.  Of  course  the  distributions  of  some  statistics  calculated 
from  samples  from  different  parents  might  be  rather  different  even 
if  the  parent  populations  are  "well"  matched  i.e.  could  not  be 
distinguished  with  available  data.  Such  statistics  should  not  be  used 
since  they  are  too  sensitive  to  parental  form.  Hence  the  athmetical 
comparisons  of  this  section  have  many  uses. 
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Kendall  (1974)  plotted  sheaves  of  von  Mises  and  Brownian 
densities  and  one  Ecose  matched  pair  of  densities.  In  the  dis¬ 
cussion  of  his  paper,  John  Kent  showed  graphs  of  Ecose  matched 
triplets  --  von  Mises,  Lack  (derived  from  a  bird  model)  and 
Brownian  densities.  The  peaks  at  6*o  and  in  the  tails  are  in 
that  order,  while  in  the  middle  of  the  range  it  reverses.  The 
degree  of  agreement  depends  upon  <  since  the  agreement  becomes 
perfect  as  k  -*■  o  and  <-►«*». 

We  wish  to  add  the  angular  Gaussian  to  those  comparisons, 
and  to  consider  the  spherical  versions  as  well.  Further  Ecose 
matching  may  be  supplemented  by  other  methods.  We  also  tried 

(a)  eye  matching  of  densities  (subjective) 

(b)  matching  so  that  the  supremum  of  the  difference  of  the 
cumulative  distributions  of  e  was  a  minimum  (objec¬ 
tive,  obtained  a  computer  search). 

We  begin  with  the  circular  versions.  As  mentioned  earlier 
Stephens  made  the  Brownian  -  von  Mises  Ecose  comparisons. 

To  match  the  circular  angular  Gaussian  to  the  von  Mises  by 
equating  values  of  Ecose,  m  and  k  must  be  related  by 


j ( sm)sds . 


From  asymptotic  expansions,  we  have 


Ij(<) 

17^7 


% 


1 


K  O, 


(26) 


K  ®  . 
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To  get  limiting  expressions  for  the  right-hand  side  of  (26)  we 
may,  for  m  -*•  o,  use 

'i2r+l 

i _  I  sm 

r(  r+2 1 r f  r+ 1  )  1 7~, 
r*o 


^(sm)  •  Z  r  |  r+2 )  r  ( r+1 )  I~  I 


and,  for  m  ■+■  ®,  use 


I.(sm)  'v  (2irsm)‘Js{l  -  3 _ }  exp  sm  . 

1  8Tm 


Some  calculations  then  yield  the  limiting  forms  of  (26) 


k  'v-  m 


<  ^  m 


df  ■ 


K  ■*  0 


(27) 


•  K 


Criterion  (a)  suggests  we  look  at  the  formulae  for  the  den¬ 
sities.  von  Hises  for  tc  -►  o  as  approximately 


(1+KCOS0) 


while  the  angular  gaussian  is,  for  m  o 


1 

In 


a 

f 


s ( l+m$cose)exp(-s  /2 )ds 


■  (1+  /7n  m  cose) , 

so  we  should  set  tc  ■  m  (ir/2)**  as  In  (27).  When  tc  Is  large, 

6  Is  trivially  shown  to  be  Gaussian  with  variance  tc"1.  When  m  -*• 
large,  tane  -  Xg/Xj  where  Xg  Is  Gaussian  (0,1)  and  X}  Is  Gaussian 
(m, 1 )  and,  since  e  is  small,  tan  e  ■  e.  Hence  e  is 
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asymptotlcal ly  Gaussian  (0,m"2).  Matching  <  ■  m2,  again  as  In 
(27).  We  In  fact  found  that  “drawing"  densities  on  a  CRT  and 
matching  the  AG(m  fixed)  and  VM  (<  to  be  found  as  tc(m))  densities 
by  eye  we  got  very  similar  ie(m)  to  the  values  got  by  Ecose  matching. 

Figure  1  shows  a  sheaf  of  angular  Gausslans.  Figure  2  shows 
a  sheaf  of  eye  matched  (Ecose  match  would  be  very  similar)  von 
Mlses  distributions.  The  relationships  may  be  seen  when  trans¬ 
parencies  of  these  figures  are  overlaid.  When  matched  pairs  are 
plotted  together  A6  exceeds  VM  at  6*0  and  x  and  Is  less  In  the 
middle  of  the  range.  With  Kent's  results,  mentioned  above,  we 
see  that  the  angular  Gaussian  is  the  extreme  member  of  the  quadru¬ 
plet  --  but  the  density  differences  are  small.  It  Is  In  fact 
hard  to  appreciate  such  pictures  (see  Figures  3b,  4b,  5b)  --  per¬ 
haps  ratios  of  densities  should  be  used. 

Matching  method  (b)  leads  us  to  more  revealing  pictures  of 
the  differences  of  the  cumulatlves  of  the  angular  Gaussian,  von 
Mlses  and  Brownian  distributions.  For  the  first  two  distribu¬ 
tions,  the  supremum  matched  k  Is,  for  all  m,  slightly  greater 
than  the  eye  -•  or  Ecose  --  matched  k  .  Of  course  the  suprema 
obtained  are  smaller  than  for  these  latter  methods.  As  is  seen 
In  Figure  3a,  the  supremum  Is  .0125  the  same  as  Stephens' value 
at  k>1.4,  when  the  Brownian  and  von  Mlses  are  compared. 

In  Figure  3b,  the  corresponding  densities  are  drawn.  Figure  4 
Illustrates  the  worst  situation  found  (at  k«2.646,  m*1.6)  for  the 
angular  Gaussian  and  the  von  Mlses  --  the  supremum  Is  less  then 
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in  Figure  3a.  In  Figure  5,  the  worst  case  for  Brownian  and 
angular  Gaussian  is  shown. 

Turning  now  to  the  sphere*  Roberts  &  Ursell  did  Ecose  match¬ 
ing  of  the  Fisher  and  Brownian  distributions  (the  matching  formula 
being  (16)  above)  and  looked  also  at  the  matching  method  (b). 

They  found  that  the  two  methods  give  very  similar  results.  Their 
conclusion  was,  of  course,  that  it  should  be  safe  to  use  the  Fisher 
distribution  for  statistics  when  the  parental  distribution  was 
really  Brownian. 

Figure  6  gives  a  sheaf  of  Fisher  densities  for  e  i.e., 
plots  of 

k  expUcose)  sine  (28) 

?sTnF7 

for  k  ■  0( . 25 ) 1 . 5 ,  2(1)10.  Figure  7  gives  a  sheaf  of  angular 
Gausslans  densities  of  e  i.e.,  plots  of 

00 

f  s2exp{-»s(s2+m2-2smcos8)  }ds  sine  (29) 

for  m  »0  (.2)1,  1.25(.25)  3.  The  matching  Is  vivid  only  when 
transparencies  of  those  figures  are  overlaid.  The  worst  case  we 
encountered  Is  shown  In  Figure  8  where  the  difference  of  the 
cumulatlves  (Figure  8a)  and  the  densities  (Figure  8b)  are  shown. 

The  supremum  In  Figure  8a  Is  slightly  less  than  .02  rather  than 
slightly  over  .02  as  Roberts  1  Ursell  found.  This  occurred  for 
a  k  of  4.56  which  corresponds  to  a  reasonably  disperse  cluster. 
Our  efforts  at  eye  fitting  of  densities  was  much  less  successful 
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here.  The  supremum  fitted  kappas  are  shown  in  the  brief  table 
below: 


m 


fitted  k 


.1 

.3 

.5 

.75 

1 

1.5 

3 

4 

5 


.160 

.482 

.815 

1.26 

1.74 

2.94 

9.4 

16.4 

25.4 


These  agree  well  with  the  formulae  obtained  from  an  Ecose 
match,  or  by  mathematically  matching  the  density  functions  (28) 
and  (29).  By  calculations  similar  to  those  used  earlier,  we  find 


k  ^  4(2ir)  _,*m  (k  *►  o) , 


K  'V* 


m 


(k  -  -)• 


(30) 


Thus  the  angular  Gaussian  Is  as  good  a  match  as  the  Brownian 


to  the  Fisher  distribution,  and  to  the  von  Rises. 
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LEGEND 


FIGURE  1.  A  sheaf  of  angular  Gaussian  densities  (formula  9) 
for  m*0( . 2 )  1.0,  1.25  (.25)  2.5  . 


FIGURE  2.  A  sheaf  of  von  Mlses  densities  (formula  1)  matched 
by  eye  to  the  curves  In  Figure  1.  Ecose  matchings  would  be 
very  similar. 


FIGURE  3.  Three  methods  for  assessing  the  best  matching  of  the 
Brownian  and  von  Mlses  distributions  when  they  are  hardest 
to  match:  v*1.06,  k«1.4  . 

(a)  A  plot  of  the  difference  BR-VM  of  the  cumulative 
distributions. 

(b)  Plots  of  the  br  and  vm  densities. 

(c)  A  plot  of  the  relative  error,  (BR-VMJ/VM  . 


FIGURE  4.  Three  methods  for  assessing  the  best  matching  of  the 
angular  Gaussian  and  the  von  Mlses  distributions  when  they 
are  hardest  to  match:  m-1.6,  k-2.646  . 

(a)  A  plot  of  the  difference  AG-VM  of  the  cumulative 
distributions. 

(b)  Plots  of  the  ag  and  vm  densities. 

(c)  A  plot  of  the  relative  error,  (AG-VM)/VM  . 


FIGURE  5.  Three  methods  for  assessing  the  best  Ecose  matching 
of  the  Brownian  and  angular  Gaussian  distributions  when 
they  are  hardest  to  match:  m-1.2,  v-0.818  . 

(a)  A  plot  of  the  difference  AG-BR  of  the  cumulative 
distributions. 

(b)  Plots  of  the  densities  ag  and  br  . 

(c)  A  plot  of  the  relative  error,  (AG-BR)/BR  . 


FIGURE  6.  A  sheaf  of  Fisher  densities  (formula  2)  for  6  for 
k»0( . 25)  1.5,  2(1)  10. 


FIGURE  7.  A  sheaf  of  angular  Gaussian  densities  (formula  8)  for 
6  for  m>0( .2)  1,  1.25  (.25)  3  . 


FIGURE  8.  Comparisons  of  the  Fisher  and  angular  Gaussian  distribu¬ 
tions  when  they  are  hardest  to  match:  ic-4.56  ,  m«2  . 

(a)  A  plot  of  the  difference  of  the  cumulatlves  for  0  ,  AG-FI 

(b)  Plots  of  the  two  densities  for  e  ,  ag3(e)  and  fl(e)  . 
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